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Abstract. We consider here the ID Grcen-Naghdi equations that are com- 
monly used in coastal oceanography to describe the propagation of large ampli- 
tude surface waves. Wc show that the solution of the Green-Naghdi equations 
can be constructed by a standard Picard iterative scheme so that there is no 
loss of regularity of the solution with respect to the initial condition. 



1. Introduction 

1.1. Presentation of the problem. The water-waves problem for an ideal liquid 
consists of describing the motion of the free surface and the evolution of the velocity 
field of a layer of perfect, incompressible, irrotational fluid under the influence of 
gravity. This motion is described by the free surface Euler equations that are 
known to be well-posed after the works of Nalimov [l6j, Yasihara [21], Craig [6], 
Wu [HO HO] and Lannes [IT]. But, because of the complexity of these equations, 
they are often replaced for pratical purposes by approximate asymptotic systems. 
The most prominent examples are the Green-Naghdi equations (GN) - which is a 
widely used model in coastal oceanography ([El El [7] and, for instance, p~8| HOj)-- 
the Shallow- Water equations, and the Boussinesq systems; their range of validity 
depends on the physical characteristics of the flow under consideration. In other 
words, they depend on certain assumptions made on the dimensionless parameters 
e, n defined as: 

a _ h% 

where a is the order of amplitude of the waves and the bottom variations; A is 
the wave-length of the waves and the bottom variations; ho is the reference depth. 
The parameter e is often called nonlincarity parameter; while (i is the shallowness 
parameter. In the shallow- water scaling (/x <C 1), and without smallness assumption 
on e one can derive the so-called Green-Naghdi equations (see [H [13] for a derivation 
and [2] for a rigorous justification) also called Serre or fully nonlinear Boussinesq 
equations [15]. 

In nondimcnsionalized variables, denoting by £(i, x) and u(t, x) the parameteriza- 
tion of the surface and the vertically averaged horizontal component of the velocity 
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at time t, and by b(x) the parameterization of the bottom, the equations read 
dtC + V • (hu) = 0, 

(h + fihT[h, eb])d t u + WC + eh(u ■ V)u 



(1) 



+ M e{ - \v[{h 3 {{u ■ V)(V ■ u) - (V ■ u) 2 )] + = 0, 



where /i = 1 + e(( — b) and 

T[h, eb]W = -jjrV(ft 3 V ' w ) + ^i^( h2 ^ b ■ w ) - h 2 \7b\7 ■ W] + e 2 V6V6 ■ W, 
while the purely topographical term ^R[h,eb]u is defined as: 

ft[h,eb]u = -L[V{h 2 {u-Vfb)-h 2 {{u-V){V ■ u) - (V ■ u) 2 )Vb] 
+e 2 ((u-W) 2 b)Wb. 

This model is often used in coastal oceanography because it takes into account the 
dispersive effects neglected by the shallow-water and it is more nonlinear than the 
Boussincsq equations. A recent rigorous justification of the GN model was given by 
Li [14] in ID and for flat bottoms, and by B. Alvarcz-Samaniego and D. Lannes [2] 
in 2008 in the general case. This latter reference relies on well-poscdness results for 
these equations given in [3] and based on general well-poscdness results for evolution 
equations using a Nash-Moscr scheme. The result of [3] covers both the case of ID 
and 2D surfaces, and allows for non flat bottoms. The reason why a Nash-Moser 
scheme is used there is because the estimates on the linearized equations exhibit 
losses of derivatives. However, in the ID case with flat bottoms, such losses do 
not occur and it is possible to construct a solution with a standard Picard iterative 
scheme as in [M]. Our goal here is to show that it is also possible to use such a 
simple scheme in the ID case with non flat bottoms, thanks to a careful analysis 
of the linearized equations. 

1.2. Organization of the paper. We start by giving some preliminary results in 
Section 12. 1\ the main theorem is then stated in Section 12.21 and proved in Section 
12.31 Finally, in Appendix [Aj we give the existence and uniqueness of a solution to 
the linear Cauchy problem associated to the Green-Naghdi equations. The proof 
of the energy conservation, stated in the main theorem, is given in Appendix [B] 

1.3. Notation. We denote by C(Ai, A 2 , ...) a constant depending on the parameters 
Ai, A2, ... and whose dependence on the Xj is always assumed to be nondecr easing. 
The notation a < b means that a < Cb, for some nonegative constant C whose 
exact expression is of no importance (in particular, it is independent of the small 
parameters involved). 

Let p be any constant with 1 < p < 00 and denote L p = L P (R) the space of all 
Lebcsgue-mcasurablc functions / with the standard norm 

l/k* = ( / l/(x)|^) 1/p <co. 



When p = 2, we denote the norm \ ■ \l 2 simply by | • I2. The inner product of any 
functions /1 and fa in the Hilbert space L 2 (R) is denoted by 

(/i,/ 2 ) = / fi(x)fa(x)dx. 
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The space L°° = L°°(R) consists of all essentially bounded, Lebesgue- measurable 
functions / with the norm 

= ess sup \f(x)\ < oo. 

We denote by W 1 ' 00 = W 1 '°°(R) = {/ E L°°,d x f € L°°} endowed with its canoni- 
cal norm. 

For any real constant s, H s = H S (M) denotes the Sobolev space of all tempered dis- 
tributions / with the norm = | vV s 1 2 < 00, where A is the pseudo-differential 
operator A = (1 - d 2 ) 1 / 2 . 

For any functions u = u(x,t) and v(x,t) defined on R x [0, T) with T > 0, we 
denote the inner product, the L p -norm and especially the L 2 -norm, as well as the 
Sobolev norm, with respect to the spatial variable x, by (u,v) = (u(-,t),v(-,t)), 
\u\lp = \u(-,t)\ L p, \u\ L 2 = \u(-,t)\ L 2 , and \u\ H * = \u(-,t)\ H *, respectively. 
Let C fe (R) denote the space of fc-times continuously diffcrentiable functions and 
C^°(R) denote the space of infinitely diffcrentiable functions, with compact sup- 
port in R; we also denote by C£°(R) the space of infinitely diffcrentiable functions 
that are bounded together with all their derivatives. 

Let / be a function of the independent variables x\, X2,---,x m ; its partial derivative 
with respect to Xk is denoted by d Xk f = f Xk for 1 < k < m. 

For any closed operator T defined on a Banach space X of functions, the commu- 
tator [T, /] is defined by [T, f]g = T(fg) - fT(g) with /, g and fg belonging to the 
domain of T. 



2. WELL-POSEDNESS OF THE GREEN-NAGHDI EQUATIONS IN ID 

For one dimensional surfaces, the Green-Naghdi equations (JTJ) can be simplified, 
after some computations, into 

f d t C + d x (hu) = 0, 

j (h + fihT[h,eb})[d t u + eud x u] + hd x ( + epihQ[h,eb]{u) = 

where h = 1 + e(C — b) and 

1 e 
T[h,eb]w = -—d x (h 3 w x ) + —[d x {h 2 b x w) - h 2 b x w x ] + e 2 b 2 x w, 

2 1 

Q[h,eb}(w) = —d x (h 3 wl) + ehw 2 x b x + e—d x {h 2 w 2 b xx ) + e 2 w 2 b xx b x . 

Remark 1. The interest of the formulation (0) of the Green-Naghdi equation is 
that all the third order derivatives of u have been factorized by (h + (j,hT[h, eb]) . 
Indeed, Q[h,eb] is a second order differential operator. This was used in [14j in the 
case of flat bottoms (b = 0). 

2.1. Preliminary results. For the sake of simplicity, we write 

1 = h + nhT[h,eb}. 
We always assume that the nonzero depth condition 
(3) 3 h > 0, inf h>h , h = 1 + e(C - b) 

is valid initially, which is a necessary condition for the GN system ([2]) to be phys- 
ically valid. We shall demonstrate that the operator T plays an important role in 
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the energy estimate and the local well-posedness of the GN system ([2]). Therefore, 
we give here some of its properties. 

The following lemma gives an important invcrtibility result on X. 

Lemma 1. Let b G C£°(R) and ( G W 1,0 °(R) be such that is satisfied. Then 
the operator 

T:H 2 (R) — >L 2 (R) 
is well defined, one-to-one and onto. 

Remark 2. Here and throughout the rest of this paper, and for the sake of sim- 
plicity, we do not try to give some optimal regularity assumption on the bottom 
parameterization b. This could easily be done, but is of no interest for our present 
purpose. Consequently, we ommit to write the dependance on b of the different 
quantities that appear in the proof. 

Proof. In order to prove the invertibility of X, let us first remark that the quantity 
\v\l defined as 

\v\l = \v\l + n\d x v\l, 

is equivalent to the iJ 1 (R)-norm but not uniformly with respect to /i G (0, 1). We 
define by i/*(R) the space 7? 1 (R) endowed with this norm. The bilinear form: 

fie . . 

a(u,v) = (hu,v) + ii{h[-j=u x —eb x u), ~^ v x ^ xV > + — (hb x u,b x v). 

is obviously continous on Hi (R) x Hi (R) . Remarking that 

a(v, v) = (Zv, v) = (hv,v) 



pie 2 



we have 



+fi(h(-^=v x - ^-eb x v) , -^=v x - -^-sb x v) + ^-(l 



i ,9 6u /. h , ,, 3e 2 ., ,, 



One deduces that 



max { 1 '^|}(l i; l2 + ' u l^^- ^Mi + ^-IMl) > \ v \l 

Since from we also get 

a(v,v) > ho\v\l + fj,h (\^v x --^-eb x v\l + ^-\b x v\i), 



h V3 

10\V\2 + yU/l0(J- 

it is easy to deduce that 



(4) a(v,v) > '^r^TW 2 



In particular, a is coercive on H%. Using Lax-Milgram lemma, for all / G L 2 (R), 
there exists unique u 6 Hl(R) such that, for all v G i?*(R) 

a(u,v) = (f,v); 
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equivalently, there is a unique variational solution to the equation 
(5) %u = f. 

Wc then get from the definition of T that 



d 



_ hu + ^d x {h 2 b x )u + e 2 nhb 2 x ii - !§d x h 3 d x u - f 



p.h 3 

3 



since u G H^R) and / G L 2 (R), we get d 2 x u G L 2 (R) and thus u G H 2 (R). □ 
The following lemma then gives some properties of the inverse operator X -1 . 

Lemma 2. Let b G C 6 °°(R), t > 1/2 and ( G iJ* 0+1 (R) be such that (0) is satis- 
fied. Then: 

(i) VO < s < t + 1, IT- 1 /!^ + v/fl^S- 1 /!*' < C(£, |/l - l\ Ht0+l )\f\„s; 

(ii) VO<s<to + l, y/]i[Z- l d x g\ B . <C(±,\h-l\ Ht0+ i)\g\ H ,; 
(Hi) Ifs>t + 1 and ( G # S (R) then: 

II ^ 1 llff=(R)^/f=(R) II ^ ||_H"-'(K)-»H S (R)< c s, 

where c s is a constant depending on \h—l\H' and independent of (fi,e) G (0, l) 2 . 

Proof. Step 1. We prove that if u G Hl(M) solves 

Tu = f + y/JLd x g 

for /, g G L 2 (R), then one has 

\u\ Hl <C(^)(\f\ 2 + \g\ 2 ). 
no 

Indeed, multiplying the equation by u and integrating by parts, one gets, with the 
notations used in the proof of lemma [1] 

a(u,u) < (f,u) - (g,^/JId x u). 

Wc thus get from the proof of Lemma [1] and Cauchy-Schwarz inequality that 



ho 



max{l, ±| } 



U|ffi < \f\2\u\ 2 + \g\2\u\ Hl , 



and the result follows easily. 

Step 2. We prove here that \%- x f\& + y/JAd^' 1 f\ H - < C(±, \h- l\ Ht0+ i)\f\ H .. 
Indeed, if / G H s and u = r £^ 1 f then %u = /. Applying A s to this identity, we get 

T(A 3 u) = A'f+[Z,A']u 
= f + VJ^d x g, 

with, 



/ = A s f - [A s , h]u + £ -^[A s ,h 2 b x }u x - e 2 fi[A s , h 2 b 



and 

~ g = YfiA s ,h i }u x ~^{A s ,h 2 b x } 



V^Us l3i„, £ V^\m 
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Now, one can deduce from the commutator estimate (see e.g Lemma 4.6 of [3]) 

(6) |[A*,F]G| 2 <|VF| Ht0 |G| ffs -i 

that 

\f\2 + \g\2 < \f\ H s+C(^-,\h-l\ H * 0+ l)(\ U \ H s-l+^jl\d x u\ H s-l). 

One can use Step 1 and a continuous induction on s to show that the inequality (i) 
holds for < s < t + 1. 

Step 3. We prove here that y /Ji\%~ 1 d x g\H= < C{-fT' \^ ~ lltf'o+Olfflff 3 ■ Indeed, if 
g £ H s and u = y /JI r Z~ 1 d x g then Tu = ^fjid x g and thus 

Z(A s u) = f + ^Jid x g, 

with, 



and 



/ = -[A s , h]u + ^[A s ,h 2 b x }u x - e 2 fi[A s , h 2 b x ]u, 



g = A s g + ^[A s ,h 3 ]u x - £ -^{A°,h 2 b x ]u. 



Proceeding now as for the Step 2, one can deduce (ii). 

Step 4. If s > to + 1 then one can prove (iii) proceeding as in Step 2 and 3 above, 
but replacing the commutator estimate ([6]) by the following one 

(7) \[A s ,F]G\ 2 <\VF\ H s-i\G\ H s-i. 

□ 

2.2. Linear analysis. In order to rewrite the GN equations ^ in a condensed 
form, let us decompose Q[h, eb](u) as 

efj,hQ[h, eb]{u) = Qi[U]u x + q 2 (U) 

where U = (£, u) T and 



(8) 



2 

Qi[U]f = -end x {h 3 u x f)+e 2 nh 2 b x u x f + e 2 fih 2 b xx uf 



q 2 (U) = e 3 ^hb xx b x u 2 + ^e 2 iid x {h 2 b xx )u 2 . 



The Green-Naghdi equations ^ can be written after applying % 1 to both sides 
of the second equation in ^ as 

(9) d t U + A[U]d x U + B{U) = 0, 

with U — (C, u) T and where 



(10) A[U] 



eu h 
T-^h-) eu + %- 1 Qi[U] 
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and 



(11) B{U) = 



eb x u 

5-^2(11) 



This subsection is devoted to the proof of energy estimates for the following 
initial value problem around some reference state U_ = (C,u) T - 



(12) 



dtU + A\U]d x U + B(U) = 0; 
U\ t = = U . 

We define now the X s spaces, which are the energy spaces for this problem. 



Definition 1. For all s > and T > 0, we denote by X s the vector space iP(R) x 
_ff s+1 (R) endowed with the norm 

VU = (C,u)eX s , \U\ 2 XS :=|C|ff. + Mk+HMff.. 

while Xf. stands for C([0, ^]; X s ) endowed with its canonical norm. 

First remark that a symmctrizcr for A[U] is given by 

/I 

(13) S=[ 

\ % 

with h = l + e{Q-b) and c l=h + fihT[h, eb) . A natural energy for the IVP H2|) is 
given by 

(14) E s (Uf = (A S U,SA S U). 

The link between E S (U) and the X s -norm is investigated in the following Lemma. 

Lemma 3. Let t e (7 6 °°(R), s > and ( G W^°°(R). Under the condition 
E S (U) is uniformly equivalent to the | • \x°-norm with respect to (//, e) € (0, l) 2 : 

E s (U)<C(\h\ L ~,\h x \ L -)\U\ X s, 

and 

\U\ X ,<C{±-)E\U). 
no 

Proof. Notice first that 

E s (Uf = \k s C\\ + (A s u,TA s u), 

one gets the first estimate using the explicit expression of %, integration by parts 
and Cauchy-Schwarz inequality. 

The other inequality can be proved by using that inf^gR h > ho > and proceeding 
as in the proof of Lemma [T] □ 

We prove now the energy estimates in the following proposition: 

Proposition 1. Let b € C£°(M.), t > 1/2, s > t + 1. Let also U = {(,u) T 
€ be such that dtU_ £ X^T 1 and satisfying the condition 0i on [0, — ]. Then for 
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all Uo £ X s there exists a unique solution U = (£, u) T £ Xj, to U2\) and for all 
< t < f 

E s (U(t)) < e eXTt E s {U ) + e f eF XT <-*- tf ')C{E a (JJ)(lf))dlf. 

Jo 

For some X T = A T (sup < t < T/e E s (U(t)), sup < t < T/e \d t h(t)\ L °°) ■ 

Proof. Existence and uniqueness of a solution to the IVP (fT!2|) is achieved in ap- 
pendix A and we thus focus our attention on the proof of the energy estimate. For 
any A £ K, we compute 

e eXt d t {e- £Xt E s (Uf) = -e\E s (Uf +d t (E s (Uf). 

Since 

E s {Uf = (A S U,SA S U), 

we have 

(15) d t (E s (Uf) = 2(A s (,A s ( t ) + 2(A s u,ZA s u t ) + (A s u, [d u ^]A s u). 
One gets using the equations (|12p and integrating by parts, 



^e eXt d t (e- £Xt E s (U) 2 ) = -y E S {U) 2 - (SA[U]A s d x U, A S U) 



-([A s ,A[U]]d x U,SA s U) - (A S B(U),SA S U) 



(16) +±(A s u,[d u %A s u). 

We now turn to bound from above the different components of the r.h.s of (JTGj) 
• Estimate of (SA[U]A s d x U, A S U). Remarking that 

eu h 



SA[U] = 



h %(su-) + Q 1 [U] 



we get 

(SA[U]A s d x U, A S U) = {euA s ( x ,A s () + {hA s u x ,A s () 

+(hA s ( x ,A s u) + ((T(eu-) + QAmWu^ A s u) 

=: A 1 +A 2 +A 3 + A A . 

We now focus to control {Aj)i<j<4. 

— Control of A\. Integrating by parts, one obtains 

A, = (euA s CA s ( x ) = -±(eu x AX,AX) 

one can conclude by Cauchy-Schwarz inequality that 

IA^KsCHuJl^E^U) 2 . 

— Control of A 2 + A3. First remark that 

\A 2 + A 3 \ = \(h x A s u,AX)\ < ll^l l^E s (U) 2 ; 



LARGE TIME EXISTENCE FOR ID GREEN-NAGHDI EQUATIONS 



9 



we get, 

\A 2 +A 3 \ <eC(\h x \ L ~)E s (U) 2 . 

— Control of A4. One computes, 

Ai = e(5(uA s u x ),A s u) + (Q 1 [U]A s u x ,A s u) 
= : An + A i2 . 

Note that 

Aa = e(h uA s u x , A s u) + ^{h 3 (uA s u x ) x , A s u x ) 

2 2 

~~Y^ b x {uA s u x ) x ,A s u) - ^(h 2 b x uA s u x ,A s u x ) 
+e 3 ^{h ub 2 x A s u x , A s u) ; 

since 

(h 3 (uA s u x ) x ,A s u x ) = ^( -(h x uA s u x ,A s u x ) + (h 3 u x A s u x ,A s u x )), 

by using successively integration by parts and the Cauchy-Schwarz inequality, one 
obtains directly: 

\An\ <£C(\u\ w i,~,\C\ w i,-)E s (U) 2 . 

For A42, remark that 

|^4 2 | = \(Qi{UWu x ,A s u)\ 

2 

= I - -efi(h 3 u x A s u x ,A s u x ) + e 2 ^(h 2 u x b x A s u x , A s u) 
+s 2 n(h 2 ub xx A s u x ,A s u)\ 

therefore 

\Ai 2 \<eC(\u\ w i,~,\(\ w i,~)E s (U) 2 . 

This shows that 

t, 

• Estimate of ( [A s , A[U}] d x U, SA S U). Remark first that 

([A s ,A[U]]d x U,SA s U) = ([A s ,eu]C x ,A s O + ([A s ,h}u x ,A s O 

+ ([A S ,T~ 1 h]Q x ,%A s u) + ([AVuK,TA s M ) 

+ ([A s ,^- 1 Q 1 [U]]u x ,U s u) 

=: B 1 +B 2 + B 3 +Bi + B 5 . 

- Control of B x + B 2 = ([A s , eu]( x , A S Q + ([A s , h]u x , A s (). Since s > t + 1, we 
can use the commutator estimate (|7|) to get 

\Bi+B 2 \ < eC(E s (U))E s (U) 2 . 



10 SAMER ISRAWI 

- Control of B4 = ([A s , eu\u Xl 2A s u). By using the explicit expression of T we get 

B 4 = ([A s ,eu}u x ,hA s u) + ^{d x [A s ,su}u x ,h 3 A s u x ) 

-y([A s , eu]u x , h 2 b x A s u x ) + ^([A s , ei A )u x , d x (h 2 b x A s uj) 

+e 2 ^i{[A s 1 eu\u x ,hb 2 x A s u), 
using the Cauchy-Schwarz inequality and the fact that 

d x [A s J]g=[A s J x ]g+[A s J]g x 

one obtains directly: 

|5 4 | < eC(E s (U))E s (U) 2 . 

- Control of B 3 = ([A s , IT 1 h\C, x ,%_A s u). Remark first that 

%[A s ,T 1 ]hQ x =%W,T 1 k]Qx - {A s ,h]( x ; 
morever, since [A S ,X _1 ] = — [A S ,2]X , one gets 

SA 8 ,^ 1 h]Q x = -[A^m^Kx + [A s ,h]( x , 
and one can check by using the explicit expression of %_ that 

SIA 8 ,^- 1 h}( x = ^^k]T 1 K x + ^d x {[A s 1 h 3 ]d x {'T 1 K x )} 

-^[A^MrtC* + y [A s , ,h 2 b x ]d x 5- 1 h<; x 

-e 2 ^[A\hb 2 x ]%- 1 hQ + [A s ,h]( x . 

One deduces directly from Lemma [21 an integration by parts, and Cauchy-Schwarz 
inequality that 

\Bs\ < C(i \h-l\ H .) {(l^l^ + f \h 2 b x \ H s+e^\hbl\ H s)\K x \ H ^ 

+(^\kl\ H *-i + £ -^\h 2 b x \ H .)\K x \H^ + \h x \ H s-AQ x \Hs-^)}\A s u\ Hi . 

Finally, since 

< C(E S (U))E S (U) and \h£\ H . + \h 2 b x \ H s < C(E S (U)), 

we deduce 

\B 3 \ <eC(E s (U))E s (U) 2 . 

- Control of B 5 = ([A s Q 1 [U]]u x ,'ZA s u) . Let us first write 

%[^,%~ l Qi\U]]^x = -[A\%%- 1 Q 1 [U]u x +[A s ,Q 1 [US\u x 
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so, that 

%[K S ,%- x QiM\n x = -[A s ,m~ 1 Qx[m^ + ^d x {[A s ,h 3 ]d x (%- 1 Q 1 \U]u x )} 

- y^A^&jT"^^} + £ -^W,h 2 b x ]d x {%- 1 Q l \U}u x ) 

-e 2 ^[A\hb 2 x ]T' 1 Q 1 [U}u x + [A s , Qi\U}W 
To control the term ([A s , Q\ [£/.]] u x , A s uj we use the explicit expression of Qi[U]'- 

Qi[U]f = \e^d(h\j) + e 2 fih 2 b x uj + s 2 fih 2 b xx uf, 
and the fact that 

d x [A s ,f]g=[A s ,d x (f-)]g. 

Similarly to control the term (d x {[A s , h ]d x (%~ Qi \U]u x )}, A s uj we use the explicit 
expression of Qi[C£], the commutator estimate ([7j and Lemma [2] Indeed, 

(d x {[A a ,h?]dx(2T 1 Qi\U]u x )}, A"u) = -|e/i([A',& 3 ]flS r (^ 1 fii e (& 3 « (II «x)),A'«x) 

-e V ( [A s , h 3 } d x {T 1 (h 2 b x u x u x ) ) , A s u x ) 
-e 2 ^([A s 1 h 3 ]d x (%-\h 2 b xx uu x )),k s u x ). 

and thus, after remarking that 

|a x (T- 1 a a! (ft 3 ti (B u a! ))| H .-i < |2 _1 ^(^^)|h» 

< II ^ 1 9 x ||j?«(E)- > H«(K) I^WffiWasliT'- 

we can proceed as for the control of -B3 to get 

\B 5 \ <eC(E s (U))E s (U) 2 . 

• Estimate of {A S B(U), SA S U). Note first that 

(eb x u 

where, q2{-) as in |(5J), so that 

(A S B(U),SA S U) = (A s (eb x u),A s O + (A s (%- 1 q 2 (U)),Z.A s u) 

= (A s (eb x u),A s C) - ( [A s , T]2T 1 g 2 © , A s u) 

+ (A s g 2 (C/),A s u ). 

Using again here the explicit expressions of 2, Q2 (LL) and Lemma [2j we get 
(A S B(|7), SA S U) < eC{E s {U_))E s {U). 
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• Estimate of (A s u, [d t ,%\A s u). We have that 

(A s u, [dt, T]A S U ) = (A s u, d t hA s u) + ^(A s u X7 d t h 3 A s u x ) 

-^{A s u,d t h 2 b x A s u x ) - ^(A s u x ,d t h 2 b x A s u) 

+e 2 ^{A s u,d t hb 2 x A s u). 

Controlling these terms by eC (E s (U_) , \dtW\L°^ )E S (U) 2 follows directly from a Cauchy- 
Schwarz inequality and an integration by parts. 

Gathering the informations provided by the above estimates and using the fact 
that H S (R) c W 1 ' 00 , we get 

e eXt d t {e~ £Xt E s {U) 2 ) < e(C(E s (U), \d t h\ L oo) - X)E S (U) 2 + eC(E s (U))E s {U). 

Taking A = A^ large enough (how large depending on sup tg j T] C(E S (U_), \ dth\L°°) 
to have the first term of the right hand side negative for all t G [0, -j], one deduces 

Vt e [0, -], e eXt d t {e~ eXt E s (U) 2 ) < eC(E s (U))E s (U). 
Integrating this differential inequality yields therefore 

Vt e [0, -], E S {U) < e eXTt E s (U ) + e [ e^^-^ClE'OJWmt' . 
e Jo 

□ 

2.3. Main result. In this subsection we prove the main result of this paper, which 
shows well-posedness of the Green-Naghdi equations over large times. 

Theorem 1. Let b G C£°(R), to > 1/2, s > to + 1. Let also the initial condition 
Uo = (C,Oi u o) T G X s , and satisfy (Q|). Then there exists a maximal T max > 0. 
uniformly bounded from below with respect to s,fj, G (0,1), such that the Green- 
Naghdi equations |lp admit a unique solution U = (£, u) T G Xj, with the initial 
condition (Co:Uq) t and preserving the nonvanishing depth condition f3J) for any 
t G [0, T "" x ) . Ln particular if T max < oo one has 

\U(t,-)\ X s — > oo as t — >T max , 

or 

inffe(t,.)=infl + e(C(*,-)-6(0) — »0 as t — ► T max . 

R R 

Morever, the following conservation of energy property holds 
d t (\t\ 2 2 + (hu,u)+Li(hTu,uj) =0, 

where T = T[h, eb\. 

Remark 3. For 2D surface waves, non flat bottoms, B. A. Samaniego and D. 
Lannes [3] proved a well-posedness result to the Green-Naghdi using a Nash-Moser 
scheme. Our result only use a standard Picard iterative and there is therefore no 
loss of regularity of the solution with respect to the initial condition. In the one- 
dimensional case and for flat bottoms, our result coincides with the one proved by 
Li in [Ip. 
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Remark 4. Our approach does not admit a straightforward generalization to the 
2D case. The main reason is that the natural energy norm X s is then given by 



which does not control the H 1 ^ 2 ) norm of u (since u takes its values in R 2 , the 
information on the rotational of u is missing). 

Remark 5. No smallness assumption on e nor fi is required in the theorem. The 
fact that T max is uniformly bounded from below with respect to these parameters 
allows us to say that if some smallness assumption is made on e, then the existence 
time becomes larger, namely of order 0(l/e). This is consistent with the existence 
time obtained for the (simpler) physical models derived under some smallness as- 
sumption on e, like the Boussinesq models. In fact, such models can be derived 
from the Green-Naghdi equations |13j . The present theorem also has some direct 
implication for the justification of variable-bottom Camassa-Holm equations [9] . 

Proof. We want to construct a sequence of approximate solution (U n = (( n ,u n )) n >o 
by the induction relation 



By Proposition [TJ we know that there is a unique solution U n+1 € C([0,oo);X s ) 
to OH) if U n e C([0,oo);X s ) and U n satisfies © for all times. Let R > be such 
that E S (U°) < R/2, it follows from Proposition Q] that U n+1 satisfies the following 
inequality 



\u\ 2 xs = \c\h + Hh+^-u\l 




E s {U n+1 {t)) < e sXTt E s {U°) + e f e eXT ^ t - t ' ) C(E s (U n {t'))dt' , 





E s {U n+1 {t)) < R/2 + (e eXTt - l)(R/2 + -f^). 



We also know from the equations that 




(18) 




Since moreover 
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we can deduce from (fT8| and the fact that h™= = 1 + e(Co — b) > ho that it is 
possible to choose T small enough for U n+1 to satisfy ([3]) on [0, — ], with ho replaced 
by ho/2. 

Finally, we deduce that the Cauchy problem 

d t U n+1 + A[U n ]d x U n+1 + B(U n ) = 0; 

\t=a " 

has a unique solution U n+1 satisfing ([3]) and the inequality 

E s {U n+1 ) < e eXTt E s (U ) + e f e eXT ^-^C(E s (U n ){t'))dt' , 



when < t < ^ and At depending only on sup te j t, E s (U n ). Thanks to this 
energy estimate, one can conclude classically (see e.g. [T]) to the existence of 

T max =T(E s (U )) >0, 

and of a unique solution U £ Xf, to @ preserving the inequality © for any 
t £ [0, T "^" ] as a limit of the iterative scheme 

rro rr , w f <9 t C/" +1 + A [£/"] £/" + 1 + B(U n ) = 0: 

U° = U , and Vn£N, J L +1 1 J " v 1 

The fact that T max is bounded from below by some T > independent of e, £ 
(0, 1) follows from the analysis above, while the behavior of the solution as t — > T ma2; 
if T maa: < oo follows from standard continuation arguments. 

Though the conservation of the energy can be found in some references (e.g. [5]), 
we reproduce it in Appendix IBl for the sake of completeness. 

□ 

Appendix A. Existence of solutions for the linearized equations 

In this section we examine existence, uniqueness, and regularity for solutions to 
the following system of equations: 

d t U + A[U]d x U = /; 
Uu. =n = U , 



(19) 



where U_ = (C,u) T G AT T is such that dtU_ £ X^ 1 and satisfy the condition ((3|) on 
[0, j]. We begin the proof by the following lemma (see for instance [17]): 

Lemma 4. Let <p £ C y °(M.) ) such that (f(r) = 1 for \r\ < 1. Let also 

J 5 = <p(5\D\), 5>0. 

Then: 

(i) Vs, s' £ R, J s : H S (R) i — ► H s ' (R) is a bounded linear operator. 

(ii) J s commutes with A s and is self-adjoint operator. 

(Hi) V/ £ C 1 (E) n L 00 ^), v £ L 2 (R) £/iere ea;iste C independent of 5 such that 

\[f,J s ]v\ Hl <C\f\cMv. 
(iv) V/ £ H S (R), s > \, J s f £ L°°(R) with 

\J S f\ L ™ < C\f\ L ~ 
where C is a constant independent of S. 
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Our strategy will be to obtain a solution to (fl9|) as a limit of solutions Us to 
(20) f d t Us + J S A[U] J 5 d x U s = /; 

For any S > 0, A 5 = J 5 A[U]J 5 is a bounded linear operator on each X s , and F = 
A s d x + B(U) E C 1 ^ 3 ) so by Cauchy-Lipschitz the ODE j20} has a unique solution, 
Us E C([0,T/e], X s ). Our task will be to obtain estimates on Us, independent of 
S E (0, 1) and to show that the solution Us has a limit as 6 \ solving (|19p . To 
do this, we remark that 

^d t E s (Us) 2 = -(SA[U]A s d x J s Us,A s J s U s ) - ([A s ,A[U}]d x J 5 Us,SA s J s U s ) 

+(A'7, SA S U 5 ) + \(A s us, [d u %A s us) 

(21) +([S, J s ]A s U s ,A[U]A s d x J s Us) + ([S, J s ]A s Us, [A s , A[U}]d x J 5 U s ) . 

Note that we do not give any details for the control of the components of the r.h.s 
(f2"Tj) other than the last two terms because the others can be handled exactly as in 
Proposition [TJ To estimate the last two terms of the r.h.s (|2"Tj) . we have that 

([S, J s ]A s Us,A[U}A s d x J s U s ) = ([l,J 5 ]A s u s ,^- 1 (hA s d x J s us)) 



-([%, J s ]A s u s ,euA s d x J l 



us) 



■ 5 us) 



+ ([% J^A'us^Qi^A'^J 
One can check by using the explicit expression of T that 

[1,J S ] - [h,J 5 ]-^d x [H,J s ]d x -^ [Hb x ,J s ]d x 

+fd x [h 2 b x ,J s }+s 2 n[hblJ s }. 

One deduces directly from Lemma |4l an integration by parts, the Cauchy-Schwarz 
inequality and the explicit expression of Qi[U] that 

([S,J 5 ]A s U s ,A[U]A s d x J s Us) < CE s (U s f, 

similarly, one can conclude 

([S, J S ]A S U S , [A s ,A\U]]d x J s U s ) < CE S (U S ) 2 , 

where C is a constant independent of 5. By the Proposition [TJ we have 

-(SA[U]A s d x J s Us,A s J s U s ) - ([A s ,A[U]]d x J 5 U s ,SA s J s Us) 

+(A'7, SA S U S ) + l(A s u Sl [d u %A s u s ) 

< CE S {U S ) 2 + CE s (f) 2 . 
Consequently, we obtain an estimate of the form 

(22) ^E S (U 5 ) 2 < CE S {U S ) 2 + CE s {ff. 

at 
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Thus Gronwall's inequality yields an estimate 

(23) E S (U S ) 2 < C(t) [E S (U ) 2 + sup E s (f) 2 } , 

[o,*] 

independent of 5 £ (0, 1). Thanks to this energy estimate, one can conclude clas- 
sically (see e.g. [T7]) to the existence of a unique solution U e C([0, T], X s ) to 

da. 



Appendix B. Conservation of the energy 

In order to prove that 

fit (ICI2 + {hu, u) + fi(hTu, uj) = 0, 

we multiply the first equation of ([2]) by £ and the second by u, integrateon R, and 
sum both equations to find 

7}dt\(\l+{d x (hu)X)+(dtu, hu)+fi(hTd t u 7 u)+(d x C, hu)+e(ud x u, hu)+ne(Qu 7 u) = 0. 
Therefore 

-fijlCl! + ^9 t (hu, u) — -(dth, u 2 ) + e(ud x u, hu) + ^(hTd t u, u) + fxe(Qu, u) = 0, 
where the term Qu is defined as: 

Qu = -^d x [(h 3 (ud 2 x u - (d x u) 2 ) + £ -[d x {h 2 ud x {ud x b) - h 2 d x b(ud 2 u - {d x uf)\ 

+e 2 hd x b{ud x (ud x b)). 
Using now the fact that h = 1 + e(( — b) and the first equation of ©. we get 

1 e 

— — (dth, u 2 ) + £(ud x u, hu) = — -(d x (hu),u) + £(ud x u, hu) = 0. 

Thus, 

(24) i5 t |C|l + ^d t (hu, u) + n(hTd t u, u) + ^(Qu, u) = 0. 

Regarding now the term /j,(hTdtU,u), we remark as in [5] that 
n(hTd t u,u) = jx(TihTid t u,u) + jx(T 2 hT 2 d t u,u), 
= \i{hT\dtu, Tiu) + fi(hT 2 d t u, T 2 u), 
= fj,(h(d t (Tiu) - d t Tiu),Tiu) + fi(hd t (T 2 u),T 2 u), 

with T*(j = 1,2) denoting the adjoint of the operators Tj given by 

T\u = —=d x u — £^—d r bu, and T 2 u = —d T bu. 
^3 2 ' - 2 

It comes: 

n(hTd t u,u) = -d t (hTiu,Tiu) - -(d t h, (Tiu) 2 ) - n(h(d t Tiu,Tiu) 

+^d t (hT 2 u,T 2 u) - !±(dth, (%uf), 

= ^dt(hT lU ,T lU ) + ^d t (hT 2 u,T 2 u) - ^(d t h, (T lU ) 2 + (T 2 u) 2 ) 
-fx(hd t Tiu,Tiu). 
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Inject this result in (j24|) to get: 

±8t(\Ql + (hu,u) + ii(hTu,u)) = ^(d t h,(T lU ) 2 + (T 2 uf) 

+fj,(hd t Tiu, Tiu) — fie(Qu, u). 
Noting that hdtTiu = dth(T\u + \/3T 2 u), it comes: 

hd t h, (T lU ) 2 + {T 2 u) 2 ) + (hd t T lU , T x u) = X - (d t h, 3(T lU ) 2 + {T 2 u) 2 + 2V3T lU T 2 uj , 

= ^(d t h,{^T lU + T 2 u) 2 ), 

where we used here the first equation of (J5J). Finally, we get: 

^d t {\C,\l + {hu,u) + n{hTu,u)^ = fie(^d x (V3T lU + T 2 u) 2 - Qu,u). 

One can easily show that (— d x (VSTiu + T 2 u) 2 — Qu, u^j = 0, which implies easily 
the result. 

Acknowledgments. The author is grateful to David Lannes for encouragement 
and many helpful discussions. 

References 

[1] S. Alinhac, P. Gerard, Operateurs pseudo-differentiels et theoreme de Nash-Moser, Savoirs 

Actucls. IntcrEditions, Paris; Editions du Centre National de la Recherche Scientifique 

(CNRS), Meudon, 1991. 190 pp. 
[2] B. Alvarez-Samaniego, D. Lannes, Large time existence for 3D water-waves and asymptotics, 

Inventiones Mathematicae 171 (2008), 485-541. 
[3] B. Alvarez-Samaniego, D. Lannes, A Nash-Moser theorem for singular evolution equations. 

Application to the Serre and Green-Naghdi equations. Indiana Univ. Math. J. 57 (2008), 

97-131. 

[4] S. V. Basenkova, N. N. Morozov, and O. P. Pogutse. Dispersiue ejects in two-dimensional 

hydrodynamics . Dokl. Akad. Nauk SSSR, 1985. 
[5] F. Chazel, Influence de la topographie sur les ondes de surface, These Universite Bordeaux I 

(2007), http// tel.archivcs-ouvcrtes.fr/tcl-00200419_v2/ 
[6] W. Craig, An existence theory for water waves and the Boussinesq and the Korteweg-de 

Vries scaling limits. Commun. Partial Differ. Equations 10, 787-1003 (1985). 
[7] A. E. Green, N. Laws, and P. M. Naghdi. On the theory of water waves. Proc. Roy. Soc. 

(London) Scr. A, 338: 43-55, 1974. 
[8] A. E. Green and P. M. Naghdi, A derivation of equations for wave propagation in water of 

variable depth, J. Fluid Mech. 78 (1976), 237-246. 
[9] S. Israwi, Variable depth KDV equations and generalizations to more nonlinear regimes. 

(2009) arXiv: 0901.3201vl. 
[10] J. W. Kim, K. J. Bai, R. C. Ertekin, and W. C. Webster. A strongly-nonlinear model for 

water waves in water of variable depth: the irrotational green-naghdi model. Journal of 

Oshore Mechanics and Arctic Engineering, Trans, of ASME,, 2003. 
[11] D. Lannes. Well-posedness of the water waves equations, J. Amer. Math. Soc. 18 (2005), 

605-654. 

[12] D. Lannes Sharp Estimates for pseudo-differential operators with symbols of limited smooth- 
ness and commutators, J. Funct. Anal. , 232 (2006), 495-539. 

[13] D. Lannes, P. Bonneton, Derivation of asymptotic two-dimensional time- dependent equations 
for surface water wave propagation, Physics of fluids 21 (2009). 

[14] Y. A. Li, A shallow-water approximation to the full water wave problem, Commun. Pure 
Appl. Math. 59 (2006), 1225-1285. 



18 



SAMER ISRAWI 



[15] Madsen, P.A. and Bingham, H.B. and Liu, H., 2002. A new Boussinesq method for fully 
nonlinear waves from shallow to deep water. J. Fluid Mech. 462, 1-30. 

[16] V. I. Nalimov, The Cauchy-Poison problem. (Russian) Dinamika Splosn. Sredy Vyp. 18 Di- 
namika Zidkost. so Svobod. Granicami,254, (1974) 104-210. 

[17] Michael E. Taylor, Partial Differential Equations II, Applied Mathematical Sciences Volume 
116. Springer. 

[18] Ge Wei, James T. Kirby, Stephan T. Grilli, and Ravishankar Subramanya. A fully nonlinear 
Boussinesq model for surface waves. w I. Highly nonlinear unsteady waves. J. Fluid Mech., 
294: 71-92, 1995. 

[19] S. Wu, Well-posedness in sobolev spaces of the full water wave problem in 2-D, Invent. Math. 

130 (1997), no. 1, 39-72. 
[20] S. Wu, Well-posedness in sobolev spaces of the full water wave problem in 3-D, J. Amer. 

Math. Soc. 12 (1999), no. 2, 445-495. 
[21] H. Yosihara, Gravity waves on the free surface of an incompressible perfect fluid of finite 

depth. Publ. Res. Inst. Math. Sci. 18 (1982), no.l, 49-96. 

Universite Bordeaux I: 1MB, 351 Cours de la Liberation, 33405 Talence Cedex, France 
E-mail address: Samer.IsrawiOmath.u-bordeauxl.fr 



